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Abstract 

A fundamental problem for the delay and backlog analysis across multi-hop paths in wireless networks is how to 
account for the random properties of the wireless channel. Since the usual statistical models for radio signals in a 

O^l propagation environment do not lend themselves easily to a description of the available service rate, the performance 
analysis of wireless networks has resorted to higher-layer abstractions, e.g., using Markov chain models. In this work, 
we propose a network calculus that can incorporate common statistical models of fading channels and obtain statistical 
bounds on delay and backlog across multiple nodes. We conduct the analysis in a transfer domain, which we refer to as 
^ the SNR domain, where the service process at a link is characterized by the instantaneous signal-to-noise ratio at the 
receiver. We discover that, in the transfer domain, the network model is governed by a dioid algebra, which we refer to 

£ — , as (min, x ) algebra. Using this algebra we derive the desired delay and backlog bounds. An application of the analysis 

C^l is demonstrated for a simple multi-hop network with Rayleigh fading channels. 

I. Introduction 

The network-layer performance analysis seeks to provide estimates on the delays experienced by traffic 
i Q | traversing the elements of a network, as well as the corresponding buffer requirements. For wireless networks, a 
question of interest is how the stochastic properties of wireless channels impact delay and backlog performance. 
Wireless channels are characterized by rapid variation of channel quality caused by the mobility and location 
of communicating devices. This is due to fading, which is the deviation in the attenuation experienced by the 
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(T) transmitted signal when traversing a wireless channel. The term fading channel is used to refer to a channel 
that experiences such effects. In this paper we explore the network- layer performance of a multi-hop network 
• where each link is represented by a fading channel. 

We model the wireless network by tandem queues with variable capacity servers, where each server expresses 
{N) the random capacity of a fading channel. We ignore the impact of coding by assuming that transmission rates 
r j~^ over the fading channels are equal to their information-theoretic capacity limit, C, which is expressed as a 
J> function of the instantaneous signal-to- noise ratio (SNR) at the receiver, 7, by C(j) = lTlog(l + 7), where 
^ W is the channel bandwidth (in Hz). Numerous models are available to describe the gain of fading channels 
!— 1 depending on the type of fading (slow or fast), and the environment (e.g., urban or rural). The instantaneous, 
information-theoretic channel capacity of a fading channel can be represented as the logarithm of 7 by (see 
Chp. 14.2 in |o)) 

C( 7 ) = clog (0(7)), (1) 

where c is a constant and the function 5(7) is used to characterize the fading channel. We are interested 
in finding bounds on the end-to-end delay and on buffer requirements for a cascade of fading channels, with 
store-and-forward processing at each channel. 

The analysis in this paper takes a system-theoretic stochastic network calculus approach [21] , which de- 
scribes the network properties using a (min, +) dioid algebra. Arrival and departure processes at a network 
element are described by bivariate stochastic processes A{t, t) and D(r,t), respectively, denoting the cu- 
mulative arrivals (departures) in the time interval [r, t). A network element is characterized by the service 
process S(r,t), denoting the available service in [r, t). The input-output relationship at the network element 
is governed by 



D(0,t) > A * S(0,t) , 



(2) 
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where the (min, +) convolution operator '*' is defined as / * g(r,t) = inf T < n <i{/(r, u) + g(u,t)}. If network 
traffic passes through a tandem of TV" network elements with service processes Si, S2, . . . , Sn, the service of 
the network as a whole can be expressed by the convolution Si * S2 * • • • * Sn- 

The stochastic properties of fading channels present a formidable challenge for a network-layer analysis since 
the service processes corresponding to the channel capacity of common fading channel models such as Rician, 
Rayleigh, or Nakagami-m, require to take a logarithm of their distributions. As discussed in the next section, 
researchers frequently turn to higherdayer abstractions to model fading channels. Widely used abstraction are 
the two-state Gilbert-Elliott model and its extensions to finite-state Markov channels (FSMC) |31|. FSMC 
models simplify the analysis to a degree that the network model becomes tractable, at least at a single node. 
Extensions to multi-hop settings encounter a rapidly growing state space. As of today, a general multihop 
analysis that is applicable to models of fading channels, such as Rician, Rayleigh, or Nakagami-m, remains 
open. 

In this paper, we pursue a novel approach to the analysis of multi-hop wireless networks. We develop a 
calculus for wireless networks that can be applied to fading channel models from the wireless communication 
literature to provide network-layer performance bounds. We view the network-layer model with arrival, 
departure and service processes as residing in a bit domain, where traffic and service is measured in bits. We 
view the fading channel models used in wireless communications as residing in an alternate domain, which we 
call the SNR domain, where channel properties are expressed in terms of the distribution of the signal-to-noise 
ratio at the receiver. We then derive a method to compute performance bounds from these traffic and service 
characterizations. 

A key observation in our work is that service elements in the SNR domain obey the laws of a dioid algebra. 
We devise a suitable dioid, referred to as (min, x) algebra, where the minimum takes the role of the standard 
addition, and the second operation is the usual multiplication, and use it for analysis in the SNR domain. In 
this domain multi-hop descriptions of fading channels become tractable. In particular, we find that a cascade 
of fading channels can be expressed in terms of a convolution in the new algebra of the constituting channels. 
The key to our analysis is that we derive performance bounds entirely in the SNR domain. Observing that 
the bit and SNR domains are linked by the exponential function, we transfer arrival and departure processes 
from the bit to the SNR domain. Then, we derive backlog and delay bounds in the transfer domain using 
the (min, x) algebra. The results are mapped back to the original bit domain to finally give us the desired 
performance bounds. Our derivations in the SNR domain require the computation of products and quotients 
of random variables. Here, we take advantage of the Mellin transform to facilitate otherwise cumbersome 
calculations. Then, the computational problem is reduced to finding the Mellin transform for service and 
traffic processes. 

The main contribution of this paper is the development of a framework for studying the impact of channel 
gain models on the network-layer performance of wireless networks. For the purposes of this paper, the SNR 
domain is used solely as a transfer domain that enables us to solve an otherwise intractable mathematical 
problem. On the other hand, the ability to map quantities that appear in network-layer models and concepts 
found in a physical-layer analysis may prove useful in a broader context, e.g., for studying cross-layer perfor- 
mance issues in wireless communications. Moreover, the (min, x) algebra and the Mellin transform form a tool 
set that can be applied more generally in wireless communications for studying the channel gain of cascades of 
fading channels. As the first paper on the (min, x) network calculus algebra, our paper only considers simple 
network scenarios and makes numerous convenient assumptions (which are made explicit in Sec. III). There 
is room for significant future work on extensions of the model and a relaxation of the presented assumptions. 

The remainder of the paper is organized as follows. In Sec. [IT] we discuss related work. We describe the 
system model in Sec. Ill, where we also motivate the use of the SNR domain. In Sec. IV we present the 
(min, x) algebra and derive performance bounds. In Sec. [V] we apply the analysis to a cascade of Rayleigh 
channels, and present numerical examples. In Sec.[Vj we investigate a network with cross traffic at each node. 
We discuss brief conclusions in Sec. IVIIL 
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II. Related Work 

Analytical approaches for network-layer performance analysis of wireless networks include queueing theory, 
effective bandwidth and, more recently, network calculus. Since the service processes corresponding to the 
channel capacity of common fading channel models such as Rician, Rayleigh, or Nakagami-m, require to 
take a logarithm of their distributions, researchers often turn to higher-layer abstractions to model fading 
channels, which lend themselves more easily to an analysis. A widely used abstraction is the two-state 



channel model developed by Gilbert 16 and Elliott |12|, and subsequent extensions to a finite-state Markov 



channel (FSMC) 135]. Markov channel models are well suited to express the time correlation of fading channel 



samples. We refer to 31 for a survey of the development and applications of FSMC models. Zorzi et al. 
|42| evaluated the accuracy of first-order Markov channel models of fading channels, where the next channel 
sample depends only on the current state of the Markov process, and higher order processes that can capture 
memory extending further back in the process history. The authors found that a first-order Markov model is a 
good approximation of the fading channel, and that using higher order Markov processes does not significantly 
improve the accuracy of the model. 

Queue-based channel (QBC) [41] is an alternative model for fading channels, which models a binary additive 
noise channel with memory based on a finite queue. Here, a queue with size M contains the last M noise 
symbols, and the noise process is an M th -oidei Markov chain. The model was found to provide a better 
approximation to the Rayleigh and Rician slow fading channels compared to the Gilbert-Elliot model [n. An 
extension of the QBC model, called Weighted QBC [l] permit queue cells (i.e., channel samples) to contribute 
with different weights to the noise process. 

Queueing theoretic studies of fading channels generally apply approximations to reduce the complexity of 
multi-hop models. Le, Nguyen, and Hossain [24] apply a decomposition approximation to analyze the loss 
probability and average delay of a multi-hop wireless network with slotted transmissions for a batch Bernoulli 
arrival process, and with independent cross traffic at each node. The wireless link is assumed to employ 
adaptive modulation and coding with multiple modes, where each mode corresponds to a given link rate, 
which is chosen based on the SNR of the channel. The channel state of a link is assumed to be stationary, and 
channel states in consecutive time slots are independent. Another decomposition approximation is presented 
by Le and Hossain |23| , who consider a multi-hop tandem network with a batch arrival process and multi-rate 
transmissions, to develop a routing scheme that can meet given delay and loss requirements. The analysis 
obtains end-to-end loss rates and delays with a decomposition analysis, and feeds the results as metrics to 
the routing algorithm. Bisnik and Abouzeid [3] model a multi-hop wireless network as an open network of 
G/G/l queueing systems. Using diffusion approximation, they obtain closed-form expressions for average end- 
to-end delays. Ishizaki and Hwang [19] studied the impact of multiuser diversity assisted packet scheduling 
on the packet delay performance in a wireless network with Nakagami-m fading channels. The network 
is modeled by a multi-queue system, where each channel is described by an FSMC. Under assumptions of 
stationarity, homogeneity, and independence of the channel processes, they approximate the tail distribution 
of packet delays, and compare it to that of a round-robin scheduler. The results indicate that the delay 
performance of multiuser diversity assisted scheduling algorithms is not necessarily superior to that of round- 
robin scheduling. Since the application of queuing theoretical methods to study cascade of fading channels 
requires many assumptions on arrival and service distributions, and simplifications of the model, the use of 
classical queueing theoretic methods for the performance analysis of multi-hop wireless networks has been put 
into question [7]. 



An effective bandwidth 22 analysis seeks to develop (asymptotic) bounds on performance metrics, e.g., 
an exponential decay of the backlog. Wu and Negi [38] have adapted an effective bandwidth analysis to the 
analysis of fading channels. They introduce the concept of effective capacity, which characterizes a wireless 
channel by a log-moment generating function (log-MGF) of the channel capacity. They obtain an asymptotic 
approximation of the delay bound violation probability of a Rayleigh fading channel. Due to the difficulty of 
computing the moment generating function (MGF) of the Rayleigh distribution, they assume non-correlated 
distributions with low SNR and estimate channel parameters from measurements. The work has been extended 



to correlated Rayleigh and correlated Nakagami-m channels, and to cascades of fading channels 36, 37 , 39 
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A closely related concept is the effective channel capacity presented by Li et al. [27], which describes the 
available channel capacity by a first order Markov chain and computes the log-MGF of the underlying Markov 
process. Taking advantage of methods developed in [26], they compute statistical delay bounds for Nakagami- 
m fading channel. Hassan, Krunz, and Matta [18] use an effective bandwidth analysis to study delay and 
loss performance at a single wireless link, which is modeled by an FSMC. For fluid On-Off traffic and FIFO 
buffering, they obtain a closed form expression for the effective bandwidth required to guarantee bounds on 
delay and packet loss. 

There is a collection of recent works that apply stochastic network calculus methods [2l] to wireless networks 
with fading channels. The stochastic network calculus is closely related to the effective bandwidth theory, 
in that it seeks to develop bounds on performance metrics under assumptions also found in the effective 
bandwidth literature. Different from effective bandwidth literature, stochastic network calculus methods seek 
to develop non-asymptotic bounds. An attractive element of a network calculus analysis is that sometimes 
it is possible to extend a single node analysis to a tandem of nodes, using the (min, +) convolution operation 
seen in the introduction. 



Fidler 14 presents a network calculus methodology for a two-state FSMC model of a single-hop fading 



channel. He applies the MGF network calculus, which was suggested in the problem sets of Chapter 7 in 
[5], and which has been developed in [13]. The MGF network calculus takes its name from the extensive 
use of moment generating functions in the derivation of performance bounds. Mahmood, Rizk, and Jiang 
|28| apply the MGF network calculus to MIMO channels and derive delay bounds for periodic traffic sources. 
Zheng et. al. 40 also use an MGF network calculus to study the performance of two-hop relay networks. A 



similar methodology is applied by Mahmod, Vehkapera, and Jiang [29] to compute the throughput of a multi- 
user DS-CDMA system with delay constraints. In the works above that apply the MGF network calculus, 
models for a cascade of multiple fading channels become complex, so that multi-node results for networks 
with more than two nodes have not been obtained. The (min, x) network calculus developed in this paper 
uses similar descriptions and assumptions for traffic and service as the MGF network calculus. By performing 
computations in a transfer domain, where fading channel models take a simpler form, we are able to compute 
multi-node service descriptions for an arbitrarily large number of nodes. 

The MGF network calculus assumes that arrivals and service at each node are independent. These as- 
sumptions can be relaxed using statistical envelope descriptions for traffic (effective envelopes) and service 
(statistical service curve) [4,21 . Jiang and Emstad [20| have applied an approach with envelopes to a fading 
channel where the wireless channel is characterized by two stochastic processes: an ideal service process and 
an impairment process, where the impairment process captures effects due to fading, noise, and cross traffic. 
Verticale and Giacomazzi [34] obtain a closed form expression for the variance of a service curve, which de- 
scribes the available service by a Markov chain. This is used for the analysis of an FSMC model of a Rayleigh 
fading channel. For computing the bounds for Markovian arrivals, they apply the bounded-variance network 
calculus introduced in [15] , which is an extension of the central limit theorem methods by Choe and Shroff 1 8 
to multi-hop paths. Verticale [33| has applied the same methodology to constant bit rate traffic. Ciucu, Pan, 
and Hohlfeld 10 and Ciucu [9] present non-asymptotic (i.e., finite number of hops), closed- form expressions 
for the delay and throughput distributions for multi-hop wireless networks. As many of the works discussed 
above, the fading channel is modeled by an abstraction that uses a link layer model of the transmission chan- 
nel. Here, the channel is assumed to be governed by a slotted-ALOHA system in half-duplex mode. The 
model of this channel is a two-state On-Off server, where a node can transmit (i.e., is in the On state) only 
when other all nodes in the interference range are not transmitting. 

There is also a literature on physical-layer performance metrics of fading channels in multi-hop wireless 
networks. Hasna and Alouini 17 have presented a framework for evaluating the end-to-end outage probability 
of a multi-hop wireless relay network with independent, non- regenerative relays, i.e., amplify-and- forward 
(AF), over Nakagami fading channels. They provide a closed-form expression for the MGF of the reciprocal 
of the equivalent end-to-end SNR for independent Nakagami fading channels. For the same AF relay network, 
Tsiftsis [32] obtains a closed-form bound for the average error probability. This bound is reportedly tight 
at low SNR, but may become loose for higher SNR values and for more severe fading environment, e.g., 
Rayleigh fading. Amarasuriya, Tellambura, and Ardakani |2j present an alternative bound to [17] on the 
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Tandem network model. 



end-to-end SNR in multi-hop AF relay networks. They derive the distribution function and the MGF for i.i.d. 
Nakagami-m fading and for independent, but non-identically distributed Rayleigh fading. The works above 
study physical layer performance bounds of channel-assisted, amplify-and-forward relaying over a multi-hop 
fading channels. They do not consider buffering or traffic burstiness, and are not concerned with network 
performance metrics addressed in this paper. Delay and backlog analysis and optimization of multihop wireless 



networks remains an open research problem 23 



III. Network Model in the Bit and SNR Domains 

We consider a wireless iV-node tandem network as shown in Fig. [TJ where each node is modeled by a server 
with an infinite buffer. We are interested in the performance experienced by a (through) flow that traverses 
the entire network and may encounter cross traffic at each node. One can think of the cross traffic at a node 
as the aggregate of all traffic traversing the node that does not belong to the through flow. The service given 
to the through flow at a node is a random process, which is governed by the instantaneous channel capacity as 
well as the cross traffic at the node. We consider a fluid-flow traffic model where the flow is infinitely divisible. 
We will work in a discrete-time domain T = {U^ : U\ = ) -U, ) € Z}, where Z is the set of integers and At is 
length of the time unit. Setting At = 1 allows us to replace tj by i, which we interpret as the index of a time 
slot. We assume that the system is started with empty queues at time t = 0. 

Different nodes and different traffic flows will be distinguished by subscripts. The cumulative arrivals to, 
the service offered by, and the departures from the node are represented by random processes A n , S n , and D n 
that will be described more precisely below, with A n = D n -\ for n = 1, . . . , N — 1. We denote by A = A\ and 
D = Dn the arrivals to and the departures from the tandem network. Throughout, we assume that arrival 
and service processes satisfy stationary bounds. 

A. Traffic and Service in the Bit Domain 

Consider for the moment a single node. Dropping subscripts, we write 

t-i t-i 
A{t, t) = ^^aj, and D{t, t) = dj , 

i=r i=r 

for the cumulative arrivals and departures, respectively, at the node in the time interval [r, t), where 
denotes the arrivals and di the departures in the i-th time slot. Due to causality, we have D(0,t) < A(0, t). 
The processes lie in the set J- of non- negative bivariate functions f(r, t) that are increasing in the second 
argument and vanish unless < r < t. The backlog at time t > is given by 

B(t) = A(0,t) - D(0,t) , (3) 

and the delay at the node is given by 

W{t) = inf {u > : A(0, t) < D(0, t + u)} . (4) 

The service of the node in the time interval [r, t) is given by a random process 5(r,t), such that Eq. ^ 
holds for every arrival process A and the corresponding departure process D. This service description with 
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bivariate functions is referred to as dynamic server. Initially defined for non-random service [6], dynamic 
servers have been extended to random processes in (5 13 



The above model is a typical network-layer model, where traffic is measured in bits, and service is measured 
in bits per second. We thus refer to this model of arrivals, departures, and service as residing in a bit domain. 
The network calculus exploits that networks which satisfy the input-output relation of Eq. ^ with equality 



can be viewed as linear systems in a (min, +) dioid algebra 25 . In the (RUj+oo}, min, +) dioid, the minimum 
and addition take the place of the standard addition and multiplication operations. The network calculus is 
based on the fact that (J 7 , min, *) is again a dioid |5|. Note that the min-plus convolution, which provides the 
second operation in the dioid, is not commutative in T . 

B. Service Model of Wireless Channel 

To compute a service model for a wireless channel, we assume that the channel state information is sampled 
at equal time intervals At. With At = 1, let ji denote the instantaneous signal-to- noise ratio observed at 
the receiver in the i-th sampling epoch. Then, 7$ is a nonnegative random variable that has the probability 
distribution of the underlying fading model. We assume that the random variables 7« are independent and 
identically distributed. This assumption is justified when At is longer than the channel coherence time. 
Otherwise, the assumption will give optimistic bounds. We emphasize that the network calculus in this 
paper applies to settings without independence, however, the derivation of performance bounds will proceed 
differently. Using Eq. 0, the instantaneous service offered by the channel in the i-th slot is given by log g("fi). 
and the corresponding service process is given by 

t-l 

S(r,t) = J2 log 9(li), (5) 

i=r 

where we haven chosen units such that the constant in Eq. ([I]) takes the value c = 1. 

The service description in Eq. §5§ requires us to work with the logarithm of fading distributions, which 
presents a non-trivial technical difficulty via the usual network calculus or queueing theory. On the other 
hand, observe that the exponential S(r,t) = e s<yT,tS> is described more simply by 

t-l 

S(r,i) = n<7( 7 • (6) 

i=T 

This motivates the development of a system model that allows us to exploit the more tractable service 
representation in Eq. ([6]). In this alternative model, arrivals, departures, and service reside in a different 
domain, where we can work directly with the distribution functions of the fading channel gain and the 
corresponding SNR at the receiver. 

C. Network Model in the SNR Domain 

We now proceed by mapping the network model from Fig. [T] into a transfer domain, which we refer to as 
SNR domain. We will seek to derive performance bounds in the transfer domain, and then map the results 
to the bit domain to obtain network-layer bounds for backlog and delays. The relationship of the network 
models in bit domain and SNR domain is illustrated in Fig. [2} 

In the previous subsection, we constructed the service process for a wireless link in the SNR domain in 
Eq. (|6j) as 

S(r,t) = e s{T ' t) . 

By analogy, we describe the arrivals and departures in the SNR domain by 

A(r,t) = e A ^, and V(t, t) = e D ^ . 

Throughout this paper, we use calligraphic upper-case letters to represent processes that characterize traffic 
or service as a function of the instantaneous SNR in the sense of Eq. ([6]). Due to the monotonicity of the 
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exponential function, T>(0, t) and -4(0, t) are increasing in t, and satisfy the causality property T>(0, t) < A(0, t). 
The backlog process is accordingly described by 

B(t) = e B(t) =A{t)/V{t). 
Since time is not affected by this transformation, the delay is given by 

W(t) = W{t) = inf{u > : A{t) < V(t + u)} . (7) 

To interpret these processes in the transfer domain, let ^ a ,i = 5 f ~ 1 (e a ' ') be the instantaneous channel SNR 
required to transmit <2j in a single time slot, assuming transmission at the rate of the capacity limit. The 
arrival process in the SNR domain can then be expressed in terms of these variables as 

t-i 

A(r,t) = ]Jg( la , i ). (8) 

i=r 

Here, we are treating channel quality expressed in terms of the instantaneous SNR as a commodity. An arrival 
in a time unit represents a workload, where j at i expresses the amount of resources that will be consumed by 
the workload. The backlog can similarly be expressed in terms of the instantaneous SNR as 

t+rs-l 

B(t)= J] 9(li), 

i=t 

with the interpretation that a node with backlog B(t) at time t requires full use of the channel capacity for 
tb time units to clear the backlog. 

Most importantly, the concept of the dynamic server translates to the SNR domain. In a network system, 
the service process in the bit domain satisfies Eq. ^ if and only if the process 5(r, t) = e s{T ^ in the SNR 
domain satisfies 



Z>(0,i)> inf {A(Q,u) -S(u,t)} . (9) 

0<u<t 

We refer to a network element that satisfies Eq. ^ for any sample path as dynamic SNR server. In this 
general setting, we not require that S takes the form in Eq. ([6]), in particular, S(r,t) does not be equal to 
S(t, u) • S(u, t). 

Traffic aggregation in the SNR domain is expressed in terms of a product. When M flows have arrivals 
at a node with arrival processes denoted by A^, k = M, then the total arrival, A^g, are given for any 
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< t < t by 

M 

A &gs {T,t)=J2 A k{T,t) ■ 
k=l 

If we let Aj and A agg denote the corresponding processes in the SNR domain, we see that 

M 

A gg (r,t) = J] A k (r,t) . 
k=i 

With the above definitions, the usual network description by a (min, +) dioid algebra in the bit domain can 
be expressed in the SNR domain by a dioid algebra on T where the second operator is a multiplication. This 



enables the development of the (min, x) network calculus in Sec. IV. We observe that the exponential function 
defines a one-to-one correspondence between arrival and departure processes in the bit and SNR domains. 
The physical arrival, departure, service, and backlog processes can be recovered from their counterparts in 
the SNR domain by taking a logarithm (see Fig. [2]). 

IV. Stochastic (min, x) Network Calculus 

This section contains our main contribution: an analytical framework for statistical end-to-end performance 
bounds for a network, where service is expressed in terms of fading distributions residing in the SNR domain. 

By an SNR process we mean a bivariate process X(t, t) taking values in R + that is increasing in the second 
argument, with X (t, t) = 1 for all t. The space of SNR processes will be denoted by J 7 " 1 ". For any pair of SNR 
processes X(t, t) and y(r, t), set 

X®y(r,t)= inf {X(r,u)-y(u,t)} , (10) 

T<U<t 

and 

M y(r,f)i„{«}. (11) 

We refer to l & and '0' as the the (min, x) convolution and (min, x) deconvolution operators, respectively. 

The arrival, departure, and service processes constructed in the previous section are SNR processes. With 
the (min, x) convolution, we can express the defining property of a dynamic SNR server from Eq. @ as 

D(0,t) > A®S(0,t) (12) 

for every pair of SNR arrival and departure processes A(r,t) and V(r,t). 

A. (min, x) Dioid Algebras 

We note that, in fact, for any system description in the bit domain by the (EU {+00}, min, +) and the 
(J-, min, *) dioid algebras there exists a corresponding characterization in the SNR domain using (Il + U 
{+00}, min, x) and (J- + , min, ®) dioids. The following argument confirm that the properties of a dioid are 
satisfied. 



Lemma 1: (R + U {+00}, min, x) is a dioid. 

Proof: We show that (R + U {+00}, min, x) satisfies the dioid axioms. For a, b, c G {R + U {+00}}: 

(1) Commutativity of mm.: min(a, b) = min(6, a). 

(2) Associativity of mm: min(min(a, 6), c) = min(a, min(fe, c)). 

(3) Idempotency of min: min(a, a) = a. 

(4) Associativity of x: (a x b) x c = a x (b x c). 

(5) Distributivity of x; min(a, b) x c = min(a x c, b x c). 

(6) Existence of a null element: The null element is +00. since min(+oo,a) = a. 
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(7) Absorption of the null element: (+00) x a = a x (+00) = +00. 

(8) Existence of a unity element: The unit of multiplication is 1, since 1 x a = a x 1 = a. 

■ 

Lemma 2: (J r+ , min, <g>) is a dioid. 

Proof: Given bivariate functions X,y,Z 6 J r+ . Since the min operation is a pointwise minimum, 
properties of the min operation, that is, properties (l)-(3) from the proof of Lemma [TJ follow from the 
(R + U {+00}, min, x) dioid. For the remaining properties we have 



(4) Associativity of ®: 



(x®y)®z{T,t) 



= inf { inf {X(r,s)-y(s,u)}-Z(u,t}} 

T<U<t T<S<U 

= inf {X(r,s)-y(s,u)-Z(u,t)} 

T<S<U<t 

= inf ' {X(t,s)- M {y(s,u) ■ Z(u,t)}} 

T<S<t S<U<t 

= x®{y®z){ T ,t) . 

(5) Distributivity of® over finite sums: 

mm(X, y ) ® Z(r,t) 

= inf {mm(X, y){r,s) ■ Z(s,t)\ 

T<S<t 

= inf {mint X{t, s) ■ Z(s,t),y(T,s) ■ Z(s,t))\ 

T<S<t 

= min (A- ®y(r,t),y <g> Z(r,t)) . 

(6) Existence of a null element: The null element is N(r,t) = +00 for all values of r and t. 

(7) Absorption of the null element: 

N®X(T,t)= inf {(+00) • X(u,t)\ = +00 . 

T<U<t 

Note that functions in F + are strictly positive by definition. 

(8) Existence of a unity element: The unity element is A(r, t), where 

A(r,t) = 

This gives 

A®X(r,t)= inf {A(T,s)-X(s,t)\ 

T<S<t 

= A(T,r)-X(T,t) 
= X(r,t). 



B. Server Concatenation and Performance Bounds 

The existing network calculus in the bit domain allows for the concatenation of tandem service elements 
using the (min, +) convolution (see page 1). As an immediate consequence, single node performance bounds 
are extended to a multi-hop setting. We now show establish the corresponding result in the (min, x ) network 
calculus. Specifically, the concatenation of dynamic SNR servers is again a dynamic SNR server. We will 
prove the result for a tandem network of two nodes, as shown in Fig. [3j 
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A(T,t) 




Z>i(r,t) = A 2 (T,t) 




Fig. 3 

Tandem of dynamic SNR servers. 



Lemma 3: Let S\(r,t) and 5 2 (r, i) be two dynamic SNR servers in tandem as shown in Fig. [3]. Then, the 
service offered by the tandem of nodes is given by the dynamic SNR server 5 ne t(r, t) with 

5 nc t(r,t) = 5i 5 2 (r,t) . 
Proof: Using Eq. the departure process V(0,t) can be written as 

V(Q,t) > inf {A 2 {0,u) -S 2 (u,t)} 

0<u<t 

> inf {inf {A(0,r)-S 1 (r,u)}-S 2 (u,t)} 

0<U<t 0<T<U 

= inf {.4(0, r) • inf {5i(r 5 u) • 5 2 (u, t)}} 

0<r<t T<u<t 

= o inf <t {4(O,r).(5 1 5 2 )(r,t)}. 

■ 

The extension to networks with more than two nodes follows by iteratively applying Lemma [3} Hence, the 
dynamic network SNR server with N nodes in tandem is given by 

5 ne t(r, t) = 5i 5 2 • • • <g> 5jv(r, i) . (13) 

Performance bounds in the (min, x ) network calculus are computed with the (min, x ) deconvolution oper- 
ator. This is analogous to role of the (min, +) deconvolution in the existing (min, +) network calculus. The 
bounds are expressed in the following lemma. 

Lemma ^: Given a system with SNR arrival process 4(r, t) and dynamic SNR server S(r,t). 

• Output Burstiness. The SNR departure process is bounded by T>(r,t) < 40 5(r, t). 

• Backlog Bound. The SNR backlog process is bounded by B(t) < 4 S(t, t). 

. Delay Bound. The delay process is bounded by W(t) < inf jd > : 4 S(t + d, t) < lj. 
Proof: For the output bound, we fix r and t with < r < t and derive 

P(0,t) 4(0, t) 

4(0, t) I 



< sup < - . 

0<u<t L .4(0, it) • S(u,t) 



rA(u,t) \ 

0<u<t ^S[U,T) J 



where we used the inequality V(0, r) > 4 5(0, r)} in the second line. 

For any fixed sample path, fix an arbitrary t > 0. The bound on the backlog is derived by 



B{t) 



4(0, t) 



— SU P 1 ^T7 ^ — 7T, — \ r 

0<n<t tA{0,U) ■ S{u,t) ) 



sup 

0<u<t 
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where we used T>(0, t) > A ® 5(0, t)} in the second step. 

Recall that the delay is invariant under the transform of domains, that is, W(t) = W(t). By definition of 
the delay in Eq. Q, a delay bound w satisfies 

W(t) = inf ( W > : . < lj 

< inf|u;>0:sup A < l\ 

~ \ ~ o< u < t \A(p,u)-S(u,t + w)J ~ / 

= inf{ W >0: sup A <lV (14) 

where we used the inequality V(0, t + w) > A 8) 5(0, t + w)} in the second line. ■ 
With an algebraic description for network performance bounds in the SNR domain in hand, we now turn 
to the problem of computing the bounds. 

C. The Mellin Transform in the SNR domain 

The concise (and familiar) expressions from the previous section for the network service and performance 
bounds in the SNR domain hide the difficulty of computing the expressions. In fact, all expressions of the 



(min, x) network calculus contain products or quotients of random variables. The Mellin transform 11 
facilitates such computations, particularly when the arrival and service processes are independent. 
The Mellin transform of a nonnegative random variable X is defined by 

M x (s) = E[X s - l \ (15) 

for any complex number s such that this expected value exists. 

Among its many properties, we will exploit that the Mellin transform of a product of two independent 
random variables X and Y equals the product of their Mellin transforms, 

Mx-y(s)=Mx(s)-M y (s). (16) 

Similarly, the Mellin transform of the quotient of independent random variables is given by 

M x /y(s) = Mx(s)-M y (2-s). (17) 

We will evaluate the Mellin transform only for sGR, where it is always well-defined (though it may take 
the value +oo). When s > 1, the Mellin transform is order-preserving, i.e., for any pair of random variables 
X,Y with Pr(X > Y) = we have Mx{s) < -My(s) for all s. When s < 1, the order is reversed. Hence 
bounds on the distribution of a random variable X generally imply bounds on its Mellin transform. 

A more subtle question is how to obtain bounds on the distribution of a random variable from its Mellin 
transform. Here, the complex inversion formula is not helpful. Instead, we will use the moment bounds 

Pr(X > a) < a~ s M x (l + s) (18) 

for all a > and s > 0. For bivariate random processes X(t, t), we will write Mx(s, r, t) = Mx(r,t)( s )- 

In our calculus, we work with the Mellin transform of (min, x) convolutions and deconvolutions, which 

not only involves products and quotients, but also requires to compute infimums and supremums. The 

computation of the exact Mellin transform for these operations is generally not feasible. We therefore resort 

to bounds, as specified in the next lemma. 

Lemma 5: Let X(t, t) and y(r,t) be two independent nonnegative bivariate random processes. For s < 1, 

the Mellin transform of the (min, x) convolution X ® y(r, t) is bounded by 

t 

M X ®y{s, T,t)<Y^ M X (s, t, u) ■ My(s, u, t) . (19) 



12 



For s > 1, the Mellin transform of the (min, x) deconvolution X y(r, t) is bounded by 

T 

Mxvy(a,T,t) < J2 M x(s,u,t) ■ My(2-s,u,r) . 



(20) 



Proof: Note that the function f(z) = : 
the convolution is estimated by 



u=0 

8-1 



is increasing for s > 1 and decreasing for s < 1. For s < 1, 



M x ®y(s,r,t) = E inf {^(r,u)-y(u,t)} 

. V T<U<t 

= E\ sup {(Xfau))'- 1 .^^))'- 1 } 

L T<1l<t. 



5-1 



In the last line, we have used the non-negativity of X and y to replace the supremum with a sum, and their 
independence to evaluate the expectation of the products. Eq. (19) follows by inserting the definition of the 
Mellin transform. The deconvolution is similarly estimated for s > 1 by 



Mxe>y(s,T,t) = E [ sup{X(u, t)/y{u, r)}) 

= e\ sup {{x{u,t)) s - x -{y{u,T)f- s } 

L 0<U<T 

T 

^Y^EKxiu^y-^-Eiiy^T)) 



s-1 



,1-si 



u=0 



and Eq. (20) follows from the definition of the Mellin transform. ■ 
As an application of Lemmas [3] and [5j we compute a bound on the Mellin transform of the service process for 
a cascade of fading channels. We also make the idealizing assumption that the channels are independent. 

Corollary 1: Consider a cascade of N independent, identically distributed fading channels, where the service 
process for each channel is given by Eq. ([6]), with i.i.d. random variables 7j. Let S ne t(r,t) denote the SNR 
service process for the entire cascade. Then, the Mellin transform for this process satisfies 



N-l + t-r 

t - T 



M Snct (s,r,t) <[" ; ' ) • (M g(l) (s)) 



t-T 



for all s < 1. 



Proof: We use the server concatenation formula in Eq. (13) to represent the service of the cascade as 



<S n et(T, t) = «->l 8) c?2 <8> ■ • • ® Sn(t, t). By Lemma [EJ its Mellin transform satisfies for s < 1 

N 



(21) 



ui ...un—i n=l 



where the sum runs over all sequences uo < u\ < • • • < un with uq = r and un = t. The assumptions on the 
service processes of the individual channels imply that each product evaluates to the same function 



N 



Y[M s (s,u n ^,u n ) = Y\{M gh) (s)) Un 



n=l 



n=l 



{M g{j) ( s )y t , 



where 7 is a random variable that has the same distribution as the 7^. Since the number of summands in 
Eq. (21) is given by (^j 1 ^ 7 ), the claim follows. ■ 
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D. Performance Bounds for the Bit Domain 

We next obtain network-level performance bounds for the bit domain. This involves a transformation from 
the SNR domain to the bit domain via the relationship in Fig. [2j which provides the translation of the abstract 
metrics T> and B into processes D and B, which, along with W, are concrete measures for traffic burstiness, 
buffer requirements, and delay. 

Theorem 1: Given a system where arrivals are described by a bivariate process A{r,t), and the available 
service is given by a dynamic server S(r,t). Let A(r,t) and S(r,t) be the corresponding SNR processes. Fix 
e > and define, for s > 0, 

min(r,t) 

U,T) . 

Then, we have the following probabilistic performance bounds. 

• Output Burstiness: Pr[D(r,t) > d £ ) < e, where 

d £ (r, t) = inf { J (log M (s, r, t) - log e) } ; 

• Backlog: Pr[B(t) > b £ ) < e, where 

b £ = inf (-(logM(s,M) -loge)); 

s>0 Is J 

• Delay: Pr{W(t) > w 6 ^ < e, where w e is the smallest number satisfying 

inf |M(s,t + w £ ,t)} < e. 

Proof: For the bound on the distribution of the output burstiness, we start from the inequality T>(r, t) < 
A® S(r,t). It follows from the moment bound and Lemma [5] that, for any choice of d > and all s > 

Pr(D(T,t) >d)= Pr{V{t) > e d ) 

< Pr(A0S{r,t) > e d ) 

< (e d )- s M A0S (l + s,T,t) 
= e- sd M(s,T,t). 

To obtain the claim, we set the right hand side equal to e, solve for d, and optimize over the value of s > 
to obtain d £ (r,t). The proof of the backlog bound proceeds in the same way, starting from the inequality 
B{t) < A(0,t)/V(0,t). 

The delay bound is slightly more subtle. Fix t > 0. Using Lemma [4] and the moment bound with a = 1, 
we obtain that 

Pr(W(t) >w)< Pr{A S(t + w,t)> 1) 

< M A 0s(s + l,t + w,t) 

for every s > 0. By Lemma [5j the Mellin transform A4a<zis{ s + l,t + w,i) satisfies a bound that agrees with 
the function M(s, t + w, t), except that the upper limit in the summation that defines M(s, t + w, t) would have 
to be replaced by r = t + w. To obtain a sharper estimate, we use instead Eq. (14) from the proof of Lemma 
|4| The resulting bound is that 

Z{t)= sup{ / (M) A 
o< u <t tS(u,t + w) y 

satisfies 

Pr(W(t) >w)< Pr(Z(t) > 1) 
<M Z (t)(s + l) 

<M(s,t + w,t), (22) 
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where we have used that the supremum in the definition of Z extends only up to u = t, and then repeated 



the proof Eq. (20). The claim follows by optimizing over s. ■ 
Corresponding bounds as in Theorem [l] can be obtained using the (min, +) algebra and the network calculus 
with moment-generating functions [13] . The significance of Theorem [T] is that it permits the application of 
the the network calculus, where traffic is characterized in the bit domain, and service is naturally expressed in 
the SNR domain. This will become evident in the next section, where the Theorem gives us concise bounds 
for delays and backlog of multi-hop networks with Rayleigh fading channels. 

V. Network Performance of Rayleigh Channels 
We now apply the techniques developed in the two previous sections to a network of Rayleigh channels. 



Consider the dynamic SNR server description for a Rayleigh fading channel, as constructed in Sec. Ill B. We 
use Eq. ([6]), with the function g{-f) given by 

5 ( 7 ) = 1 + 7 = 1 + 7H 2 , (23) 

where 7 is the average SNR of the channel and \h\ is the fading gain. For Rayleigh fading, \h\ is a Rayleigh 
random variable with probability density f(x) = 2xe~ x . In a physical system, 7 = P r /a 2 , where P r and 
a 2 are the received signal power and the (additive white Gaussian) noise power at the receiver, respectively. 
Then, \h\ 2 is exponentially distributed, and the Mellin transform of 5(7) is given by 

M gh) (s) = e 1 /^ S - 1 T( S ,T 1 ), 



where T(s, y) = x s 1 e x dx is the upper incomplete Gamma function. Using the assumption that the 7 
are independent and identically distributed, we obtain for the Mellin transform of the dynamic server 



M s (s,r,t) = (e^V^r^r 1 ))" '• (24) 



For the arrival process, we use a characterization due to Chang [5J, where the moment-generating function 
of the cumulative arrival process in the bit domain is bounded by 



- log E[e sA ^] < p(s) ■ (t - t) + a(s) 



for some s > 0. In general, p(s) and a(s) are nonnegative increasing functions of s that may become infinite 
when s is large. This traffic class, referred to as (cr(s), p(s)) bounded arrivals, is broad enough to include 
Markov-modulated arrival processes. The corresponding class of SNR arrival processes is defined by the 
condition that 

M A (s, t, t) < e (s-iy( P (s-iy(t- T )+a(s-i)) ( 25 ) 

for some s > 1. 

A . Performance Bounds of Rayleigh Fading Channels 

We consider the transmission of (a(s), p(s)) bounded arrivals on a Rayleigh fading channel. To obtain 
single-hop performance bounds, we apply Theorem [T] with the expressions for the Mellin transforms for the 
SNR service and arrival processes from Eqs. (24) and ( |25[ ). For the function M(s,r, t) from the statement of 



the theorem, we compute 

00 

M(a,T,t) < e <P(s)(t-T)+<r(s)) V- (e»"We 1 /^-«r(l-a > 7- 1 )) tt , 



u=[r-t]+ 



V(s) 
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where [r — 1]+ is the maximum of r — t and 0. The sum converges when V(s) < 1, which can be interpreted as 
a stability condition. Inserting the result into Theorem [TJ we obtain for the output burstiness the probabilistic 
bound 

d £ net (r,t) = mf{p( S )(t - r) + a(s) - ~(log(l- V( S )) + loge)} . 
The backlog bound is obtained by setting r = t, 

b £ = inf lots) - -(logfl - V(a)) + loge) 1 . 

s>0 { S J 

The delay bound is the smallest number w £ such that 

5 f .{ i-v ( .) 

We also derive end-to-end bounds for a cascade of N Rayleigh channels with (cr(s), p(s)) bounded arrivals, 
using the same parameters as before. Let S ne t(r, t) be the service process for the entire cascade. By Corollary[TJ 
its Mellin transform satisfies for < r < t 

*!*.(., r,l) < ( N ~^*~ T ) ■ (eW-TTCT 1 ))^. 
We will use again Theorem [Tj For < r < i, we compute 

e s-(p(s)(t-r)+a(s)) 

M net (s,r,t)< {1 _ v(s))N , 
where V(s) is as defined above, and where we applied the combinatorial identity 

u=0 v / \ ) 

for any iV > 1 and for < x < 1. Inserting M ne t(s,T, t) into Theorem [l] gives for the end-to-end output 
bound, denoted by d^ et (r, t), the value 

< et (r,i) = inf " r) + - ^(iVlog(l-F( S )) + loge)} . 

Note that for N = 1, this bound agrees with the previous bound for a single node. In the same way, we derive 
the probabilistic backlog bound 

b £ net = inf - J(iVlog(l - V(s)) + log E )| . (27) 

For the delay bound, we estimate for w > 

M net (s, r + t) < e-^PW^-W) ^ f N ~ 1 + ^ (V{s)) u 

u=w ^ ' 

s(-p(s)«)+o-(s)) 

- (1-^,))" • (28) 

Here, the first term in the minimum is obtained by extending the summation down to u = 0, and the second 
term results from the inequality 

N -l + u\ at-, [N-l + u-w 
< w 



U J \ u — w 
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for u > w. In both cases, the resulting sum can be evaluated with Eq. (26). The delay bound w £ is determined 
according to Theorem[T]by setting the right hand side of Eq. (28 ) equal to e, solving for w, and minimizing over 
s. Because of the complexity of the bound in Eq. (28), the last two steps can only be performed numerically. 

It is apparent that the complexity of computing end-to-end bounds is no different than bounds for a 
single channel. More importantly, observe that the end-to-end bounds scale linearly in the number of nodes 
N. Relaxing the strong independence assumptions on the channel properties would give different scaling 
properties. 

B. Numerical Examples 

We next present numerical results, where we assume Rayleigh channels with a transmission bandwidth 
of W = 20 kHz. For traffic, we use (<r(s),p(s)) bounded arrivals with default values o~(s) = 50 kb and 
p(s) = 30 kbps, i.e., the bounds on rate and bursts are deterministic. Hence, the source of randomness in the 
examples results only from randomness of the channels. We use a violation probability of e = 10~ 4 . 

In Fig. [4] we show the end-to-end backlog for cascades of N Rayleigh channels, as a function of the average 
SNR of each channel. Even though the backlog bounds increase only linearly in the number of nodes, the 
per-node requirements - at least for the last node of the cascade - must satisfy the end-to-end bounds, since 
it cannot be assumed that backlog is equally distributed across the nodes. When the SNR of the nodes is 
sufficiently high, the backlog remains low even for a large number of hops. We observe that the channel 
becomes saturated for 7 = 5 dB. When the number of channels is small, the backlog increases sharply in the 
vicinity of 7 = 5, but remains low everywhere else. 

In Figs. [5] and [6] we present, for values of the average SNR of 7 = 10 dB (Fig. [5]) and 7 = 20 dB (Fig. [6]), 
how the end-to-end backlog increases as a function of the transmission rate, for different network sizes. We 
observe that the maximum achievable rate that does not result in a 'blow-up' of the backlog decreases as the 
number of nodes is increased. We observe that the blow-up occurs earlier when the channel has a lower SNR. 

Suppose that buffer sizes are set to satisfy the end-to-end backlog. For a fixed buffer size 6 max , we can then 
use the probability Pr(B nct (t) > 6 max ) as an estimate of the probability of dropped traffic, and refer to it 
as the loss probability. In Figs. [7] and [8j we depict the loss probability as a function of the average channel 
SNR, for values of 6 max = 200 kb (Fig. Q and 6 max = 400 kb (Fig. g), for traffic with a rate of p(s) = 20 and 
30 kbps, and for N = 1, 10, and 20 nodes. The figure shows the minimum SNR needed to support a given 
loss probability is very sensitive to the number of network nodes. 

In Figs. [9] and [TU] we show the violation probability for given end-to-end delay bounds, where we compare the 
delays at a single node (N = 1) with a multi-hop network (N = 10) for different average channel SNR values, 



where we use Eqs. (22) and (|28[). The traffic parameters are o~(s) = 50 kb for the burst, and p(s) = 20 kbps 
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Fig. 5 

End-to-end backlog bound {b e net ) vs. arrival rate (p(s)) for multi-hop Rayleigh fading channels with 
e = 10~ 4 , (a(s),p(s)) bounded traffic with <t(s) = 50 kb and 7 = 10 dB, and W — 20 kHz. 




(in Fig. [9]) p(s) = 30 kbps (in Fig. 10). The figures illustrates that at sufficiently high SNR values, low delays 
are achieved even when traffic traverses 10 links. When the SNR is decreased, we can observe how the delay 
performance deteriorates in the multi-hop scenario. Note that the curves for the violation probability are 
essentially piecewise linear with two segments. This is caused by the different exponential decay rates, which 
follows from Eq. (28). Depending on where in the equation the minimum occurs, we obtain a faster or slower 
decay. 

Next we investigate how the delay violation probability for given end-to-end delay bounds is impacted by 
the SNR of the channel. We consider a network with N = 10 nodes, with the same parameters as before. 
The traffic parameters are a(s) = 50 kb and p(s) = 20 kbps. We consider end-to-end delay bounds of 
w £ = 10,20,30,40 ms. Figs. \TT\ presents the results. An interesting observation is that the SNR required to 
support a given violation probability for a short delay bound w £ = 10 ms, the SNR requirement increase fast 
for low violation probabilities. Delay bounds w e = 20 ms and higher can be supported with low violation 
probabilities even when the average SNR is small. 

The graphs presented here can be used in the planning of a multi-hop wireless network where predefined 
QoS bounds are desired for a given transmission rate. Since the average SNR depends largely on the path loss, 
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Fig. 7 

LOSS PROBABILITY (s(b)) VS. AVERAGE CHANNEL SNR (7) FOR MULTI-HOP RAYLEIGH FADING CHANNELS FOR 
TV = 1,10 AND 20, WITH BUFFER SIZE 200 KB, (a(s),p(s)) BOUNDED TRAFFIC WITH a(s) = 50 KB AND p(s) = 20 OR 

30 kbps, and W = 20 kHz. 
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Fig. 8 

LOSS PROBABILITY (e(b)) VS. AVERAGE CHANNEL SNR (7) FOR MULTI-HOP RAYLEIGH FADING CHANNELS FOR 
N = 1,10 AND 20 WITH BUFFER SIZE 400 KB, (cr(s), p(s)) BOUNDED TRAFFIC WITH a(s) = 50 KB AND p(s) — 20 OR 

30 kbps, and W = 20 kHz. 



which, in turn, is a function of the transmission radius, the graphs could help with determining the maximum 
distance between nodes to support a desired transmission rate and QoS. 
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Fig. 9 

Delay bound violation probability (e(w)) vs. end-to-end delays for multi-hop Rayleigh fading 
CHANNELS FOR N = 1, 10, 7 = 5, 10, 15,20 dB, (a(s),p(s)) bounded traffic with er(s) = 50 KB AND 

p(s) = 20 kbps and W = 20 kHz. 
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Fig. 10 

Delay bound violation probability (e(w)) vs. end-to-end delay bound for multi-hop Rayleigh fading 
channels for N = 1, 10, 7 = 5, 10, 15, 20 dB, (a(s),p(s)) bounded traffic with a(s) — 50 kb and p(s) = 30 

kbps and W = 20 kHz. 
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Fig. 11 

Delay bound violation probability (e(w)) vs. average channel SNR (7) for multi-hop Rayleigh 

FADING CHANNELS FOR TV = 10, W £ = 10,20,30,40 MS, (a(s),p(s)) BOUNDED TRAFFIC WITH cr(s) = 50 KB AND 

p(s) = 20 kbps and W = 20 kHz. 
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Fig. 12 

Single Fading Channel with Cross Traffic. 



VI. Fading Channels With Cross Traffic 



Consider a scenario in Fig. 12 where a through flow arriving to a fading channel shares the available 
bandwidth with other flows. We will refer to the traffic from these other flows as cross traffic. We use A d (t, t) 
and A c (t, t) to denote the SNR arrival processes of the through flow and the cross traffic, respectively, and 
let T> (t, t) and V c (r,t) denote the corresponding departure processes. In the SNR domain, cross traffic can 
be viewed as reducing the channel capacity of the through flow by generating interference. 

The following lemma states that, in the SNR domain, the service available to a through flow that experiences 
cross traffic at a channel can be expressed by a dynamic SNR server. 



Lemma 6: Consider a channel with a through flow and cross traffic as shown in Fig. 12 Assume that the 
channel provides a dynamic SNR server to the aggregate of through flow and cross traffic, with service process 
S(r,t), i.e, 

V o (0, t) ■ D c (0, t) > {A ■ A c ) ® 5(0, t) . 

Then 

S ^ t) = A^J) 
is a dynamic SNR server satisfying for all f > that 

X>o(0,t) >Ao®S o (0,t). (29) 

If, moreover, the service to the cross flows satisfies the upper bound T> c (0, t) < A c <8> 5(0, t), then 

T> o (0,t) > A ® max{l,5 o }(0,t) . (30) 
Proof: For any a sample path, and any t > 0, we have 

V o {0, t) ■ P c (0, t) > inf {(A(0, r) • A c (0, r)) • 5(r, t)} . 

0<T<t 

Let t* be the point where the infimum is assumed. Dividing by T> c (r,t), we obtain 

(Ao • Ac) ®S(0,t) 



V o (0,t) > 



> 



v c (o,t) 

A o (0,T*)-Ac(0,T*)-S(T*,t) 



v c (o,t) 



where we used that T> c (0, t) < A c (0,t) by causality. This the first claim in Eq. (29). To see the second claim, 
assume that £> c (0, t) < A c <8> 5(0, t). Then V c (0, t) < A c (0, t*) ■ S(r*,t), and therefore 

A o (0,T*)-A c (0,T*)-S{r*,t) 



V o (0,t) > 



A c (p,r*)-S(r*,t) 

A {T*,t). 
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t-T 



Combining this with the first part of the proof, we obtain 

Z>o(0,t) > A(0,r*) • max{l,S (r*,t)}, 

proving the second claim. ■ 
Note that S (r,t) need not be monotone in i, i.e., it may not lie in T + . 

A. Performance Bounds With Cross Traffic 

We next estimate the service process available to the through flow across a cascade of channels. Assuming 
that A c and S are independent, we obtain for the Mellin transform of the service process at a single node 

M So (s,T,t) = M s/Ac (s,T,t) 

= M s (s,T,t)-M Ac (2-s,T,t). (31) 

These service descriptions may be convolved, using Lemmas [3] and [5] to obtain a bound for the Mellin transform 
of the service provided by a cascade of fading channels to a flow that experiences cross traffic at each node. 

Corollary 2: Consider a network as in Fig. [TJ where a through flow experiences cross traffic at each fading 
channel. Let the SNR service process at each channel be given by Eq. ([6]). Assume that the arrival process of 
the cross traffic satisfies Eq. (25) with parameters (<7 c (s), p c (s)). Assume that arrivals from through flow and 
cross traffic, as well as the service processes at each channel are independent. Then, the service provided by 
the network to the through flow satisfies V o (0, t) > A ® 5 o ,net(0, t), where the Mellin transform of 5 0jne t(r, t) 
satisfies, for s < 1, 

M So , a Js,r,t) < e (i-)-^ B (i-.)^-ljH-^ . (AWs ) e a-s)-Mi-s)) 
Proof: For a single node (N = 1), we estimate for < r < t and s < 1 
M So (s, r, t) = M s {s, r, t) ■ M Ac (2 - s, r, t) 

where we have used that 2 — s > 1 for s < 1 . By Lemma [3j the service of the network is given by the 
(min, x ) convolution, S 0>ne t(T, t) = <S 0j i ® . . . ® 5 0) at(t, t). We use Lemma [5] to bound its Mellin transform by 

N 

M Snct (s,T,t) < Yl M So, n ( s ^n-l,U n ) , 

U\...UN -1 71=1 

where the sum runs over all sequences uq < U\ < ■ ■ ■ < un with uq = r and un = t. Collecting terms, we 
compute, as in the proof of Corollary [TJ that each product evaluates to the same term, 

N 

II M So Js,u n . u u n ) = e (i-)^(H(^ w ( s) . e d-)-0%(i-.))*- . 

n=l 

Since the number of summands is ( N ~ < 1 J t T < ~ T ) , this proves the claim. ■ 
We now consider Rayleigh fading channels, and assume that both through flow and cross traffic are 
(a(s), p(s)) bounded traffic with parameters a Q (s) and p (s) for the through flow, and a c (s) and p c (s) for 
the cross traffic. We next give end-to-end performance bounds, using Theorem [TJ For the function M(s,t, t), 
we compute as in Sec. [V}B for < r < t, 



e s-(p (s)(t-T)+a {s)+Na c (s)) 

where 



M 0i net(s,T,t) < , , , • ; v 



V (s) = e s '(^( s )+^W e 1 /7^-^ r ( 1 _ 
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This computation is valid under the stability condition that V (s) < 1. Thus, we obtain that the output 
burstiness at the network egress gives 



d £ ,net(r,t) < inf I a a {s) + Na c {s) + Po {s)(t - r) - - 

s>0 s 



Nlag(l-V (s)) +loge 



For the end-to-end backlog of the through flow, we obtain 



b £ ,nct(t) < mf\a (s)+Na c (s)-- 

s>0 S 



N log (I -V (s)) +loge 



Similarly, for the delay bound, we estimate for w > 



M net ( S ,t + W,t) < e s(-Po(s) W +a (s)+Na c (s)) J- 1 + 

u=w ^ 



{ p s(-p (s)w+cr (s)+Ncr c (s)) 
(i-v.(.))» -('•«W» V "1 



B. Numerical Examples 

We now present numerical examples for Rayleigh fading channels. We use the same traffic and channel model 
as in the numerical examples in Sec. V-B . The cross traffic is (a(s), p(s)) bounded traffic with parameters a c (s) 
and p c (s). As before, we assume that cross traffic is deterministically bounded by using fixed values for a c (s) 
and p c {s). Consequently, there is no statistical multiplexing between through and cross traffic. Throughout, 
the violation probability is set to e = fCP 4 . 

In Fig. 13 we show the end-to-end backlog bound b £ D net as a function of the average channel SNR 7. The 
ow has fixed parameters a {s) = 50 kb and p (s) = 30 kbps, and the cross traffic has parameters 
50 kb and p c {s) = 0, 10,20. We consider networks with iV = 10 and iV = 100 nodes. The graph 



through fC 

illustrates how the bursts of the cross traffic contribute to the backlog bound. There is an additive component 
for each traversed channel, which explains the difference for the backlog with 10 and 100 channels. The 
minimum required SNR needed for stability appears less sensitive to the number of channels. 

In Figs. 14 and 15 we again evaluate the end-to-end backlog bound b £ onet . Here, we keep the channel SNR 
values constant at 7 = 10 (Fig. 14) and 7 = 20 (Fig. 15). We vary the rate of the through flow p {s) on the 
x axis, and plot graphs for different values of p c (s) = 0, 10, 20. We again consider iV = 10 and iV = 100. The 
outcomes are as seen in the previous figure. We observe the effects of the additive component contributed by 
each node. Also, we observe, at least for lower transmission rates, that the stability of end-to-end backlogs is 
not sensitive to the number of traversed channels. 
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Average channel SNR (dB) 



Fig. 13 

End-to-end backlog bound (b e o nct ) vs. average channel SNR (7) for multi-hop Rayleigh fading 

CHANNELS WITH £ — 1CT 4 , N — 10 AND 100, (a(s),p(s)) BOUNDED TRAFFIC WITH CT (s) = 50 KB, p {s) — 30 KBPS, 
(a(s),p(s)) BOUNDED CROSS TRAFFIC WITH (T c (s) = 50 KB AND p c (s) — 0,10,20, AND W = 20 KHz. 
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10 nodes, p =0 
10 nodes, p =10 
10 nodes, p =20 
100 nodes, p = 
100 nodes, p = 10 
100 nodes, p =20 



30 40 
Arrival rate (kbps) 



Fig. 14 

End-to-end backlog bound (b e onct ) vs. arrival rate (p {s)) for multi-hop Rayleigh fading channels 
with e = 10~ 4 , (a(s),p(s)) bounded traffic with a (s) — 50 KB, (<j(s),p(s)) BOUNDED cross traffic with 
cr c (s) = 50 kb and p c (s) = 0, 10, 20, 7 = 10 dB, and W = 20 kHz. 




20 



40 60 80 

Arrival rate (kbps) 



100 



120 



Fig. 15 

End-to-end backlog bound (b e o nct ) vs. arrival rate (p (s)) for multi-hop Rayleigh fading channels 
WITH £ — 10 -4 , (a(s),p(s)) BOUNDED traffic with <J (s) = 50 KB, (a(s) 1 p(s)) BOUNDED cross traffic with 
cr c (s) = 50 kb and p c (s) = 0, 10, 20, 7 = 20 dB, and W = 20 kHz. 
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VII. Conclusion 

We have developed a novel network calculus that can incorporate fading channel distributions, without 
the need for secondary models, such as FSMC. We use the calculus to compute statistical bounds on delay 
and backlog of multi-hop wireless networks with fading channels. We took a fresh point of view, where the 
descriptions of the arrivals and the fading channels reside in different domains, referred to as bit domain and 
SNR domain. We found that by mapping arrival processes to the SNR domain, an end-to-end analysis with 
fading channels becomes tractable. An important discovery was that arrivals and service in the SNR domain 
obey the laws of a (min, x) dioid algebra. 

The analytical framework developed in this paper appears suitable to study a broad class of fading channels 
and their impact on the network-layer performance in wireless networks. Even though we computed numerical 
examples for very simplified networks, in particular, we made strong independence assumptions for the fading 
channels, our (min, x) network calculus is applicable to networks where these assumptions are relaxed. Gen- 
eralizing our framework and obtaining a more profound understanding of the dioid algebra and computational 
methods in the SNR domain is the subject of future research. 
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